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> ' Abstract 

00 ' 

pg ■ We compute the gluon distribution in deep inelastic scattering at small x by solv- 

ing numerically the angular ordering evolution equation. The leading order con- 
tribution, obtained by neglecting angular ordering, satisfies the BFKL equation. 
Our aim is the analysis of the subleading corrections. Although not complete - 
Q_i! the exact next-to-leading contribution is not yet available — these corrections are 

important since they come from the physical property of coherence of QCD radia- 
tion. In particular we discuss the subleading correction to the BFKL characteristic 
function and the gluon distribution's dependence on the maximum available angle. 
Conformal invariance of the BFKL equation is lost, however this is not enough to 
bring the small- x gluon distribution into the perturbative regime: although large 
momentum regions are enhanced by angular ordering, the small momentum regions 
are not fully suppressed. As a consequence, the gluon anomalous dimension is finite 
and tends to the BFKL value 7 = 1/2 for as — > 0. The main physical differences 
with respect to the BFKL case are that angular ordering leads to 1) a larger gluon 
anomalous dimension, 2) less singular behaviour for x — > and 3) reduced diffusion 
in transverse momentum. 
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1 Introduction 



Angular ordering is an important feature of perturbative QCD [T|] with a deep theoretical 
origin and many phenomenological consequences ||. It is the result of destructive interference: 
outside angular ordered regions amplitudes involving soft gluons cancel. This property is quite 
general. It is present in both time-like processes, such as e + e~ annihilation, and in space-like 
processes, such as deep inelastic scattering (DIS). Moreover it is valid in the regions both 
of large and small x, in which x is the registered energy fraction in the e + e~ fragmentation 
function or the Bjorken variable in the DIS structure function. Due to the universality of 
angular ordering one has a unified leading order description of all hard processes involving 
coherent soft gluon emission. 

Angular ordering is important in the calculation of multi-parton distributions by resum- 
mation of powers of InQ 2 , with Q the hard scale, and of powers of In a; or ln(l — x) for small 
or large x. This is due to the fact that angular ordering defines the structure of the collinear 
singularities and, to leading order, their relation to the infra-red (IR) singularities for x — > or 
x — > 1. In particular one finds that collinear singularities in the emitted transverse momenta 
contribute both to InQ and In a; or ln(l — x). This is because angular ordering implies ordering 
in the emitted transverse momenta divided by the energies. 

In this paper we start a systematic study of multi-parton emission in DIS at small x. The 
detailed analysis of angular ordering in multi-parton emission at small x and the related virtual 
corrections has been done in Ref. M (see also ||), where it was shown that to leading order the 
initial-state gluon emission can be formulated as a branching process in which both angular 
ordering and virtual corrections are taken into account. In this first paper we study the fully 
inclusive gluon density which gives the structure function at small x. This gluon density is 
given by an inclusive recurrence equation deduced from the small-a; coherent branching (the 
CCFM equation). 

In spite of the universality of angular ordering, the space-like and time-like distributions in 
the small- a; region are profoundly different even to leading order. In e + e~ annihilation processes 
the small-a; distributions are obtained by resumming the In a; powers which come both from 
IR and from collinear singularities in the angular ordered regions. In DIS, angular ordering is 
essential for describing the structure of the final state, but not for the gluon density at small-a;. 
This is because in the resummation of singular terms of the gluon density, there is a cancellation 
between the real and virtual contributions. The only remaining collinear singularity is the one 
originating from the first gluon emission. As a result, to leading order the small-a; gluon density 
is obtained by resumming In x powers coming only from IR singularities, and angular ordering 
contributes only to subleading corrections. 

The calculation of the gluon density by resummation of In x powers without angular ordering 
was done 20 years ago |5[| and leads to the BFKL equation for J-(x, k), the gluon density at fixed 
transverse momentum k, related to the small-a; part of the gluon structure function F(x, Q) by 

F{x,Q)= fdk 2 F(x,k)Q{Q-k). (1) 



In the moment representation the BFKL equation has solutions of the form 

du /l\ w _ ,., _ ,., 1 (k 2 ^ 1 



Xj7{x > k) = ' M\x) * u{k) ' ^(*)~^(f) ( 2 ) 



with fco an arbitrary constant and 7 given by a solution of the well known BFKL characteristic 
equation 

1 = — X(7) , X(7) = 2^(1) - ^(7) - ^(1 - 7) , (3) 

where 0:5 = A " s and ^ is the logarithmic derivative of the gamma function. The QCD 
coupling as is taken as a fixed parameter. The renormalisation group dependence of as on a 
scale is an effect which goes beyond this leading order contribution in which one resums the 
powers (as/uj) n for u — > 0. The next-to-leading order contribution, as(as/u>) n , is so far only 
partially computed ||. 

In eq. ([!]) 7 plays a role analogous to that of the gluon anomalous dimension, however its 
origin is not the renormalisation group but conformal invariance, which is a consequence of the 
absence of collinear singularities and which implies that in !F(x, k) all regions of k 2 large or 
small are equally important. This is reflected in the symmetry 7 — > 1 — 7 of the characteristic 
function. 

The gluon distribution in the angular ordered equation depends on an additional variable 
p, which corresponds to the maximum available angle for the initial state radiation. We denote 
by Au(k,p) the gluon distribution in this case. The solutions have a form similar to (0) 



^(M~^(|j G(p/k), (4) 

where G(p/k), describing the angular dependence, has a structure typical of a form factor, 
vanishing for p —>■ 0. As in (|3p, the gluon anomalous dimension 7 is given by a modified char- 
acteristic function xili a s) which depends also on as- Angular ordering breaks the conformal 
invariance of the gluon density, so that the modified characteristic function is not symmetric 

for 7 — » 1 — 7. 

To leading order, i.e. for the leading powers ag/u n , the two equations are equivalent, so that 
*Ao(k,p) — > Tu}{k). We have therefore that the angular dependence in G{p/k) is a subleading 
correction of order as- The gluon anomalous dimension has the expansion 



7 = 7 + a 5 7i H (5) 

where the leading order result 7, given by the BFKL characteristic function, and the first 
correction 71 are functions of the ratio as/oo. We have that as'ji is the part of the next-to- 
leading correction of the gluon anomalous dimension which comes from angular ordering. 

The study of the differences between the two equations will be done by analytical and 
especially by numerical calculations. As we shall discuss, the equation with angular ordering 
can be diagonalised only partially thus numerical methods are needed. The study of some of 
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Figure 1: Kinematic diagram for a DIS process at parton level: k- m is the incoming gluon, 
defined to have energy E; the fcj are the exchanged gluons (k n is the gluon which undergoes 
the hard collision) and the qi are the gluons emitted in the initial state. 



the phenomenological features of the gluon density with angular ordering has been done in Ref. 

0- 

In Sect. 2 we recall the main elements of the CCFM equation and its relation to the BFKL 
one. We discuss how a hard scale enters. We deduce some simple analytical properties and 
the behaviour of the solution. In Sect. 3 we discuss the numerical methods used. In Sect. 4 we 
present the results. Sect. 5 contains a summary of the main points and some conclusions. 



2 Equation for the gluon density 

In this section we recall the basic elements for the small-x coherent branching and the inclusive 
equation for the gluon density. 

We start by considering the kinematic diagram for a DIS process at parton level represented 
in fig. 1. All partons involved are gluons since gluons dominate the small- x region. The last 
exchanged parton n undergoes the hard collision at the scale Q. For the exchanged gluon i 
we denote by Xi and fcj the energy fraction and transverse momentum with respect to the 
incoming gluon. Introducing the energy ratio Z{ = x^/x^i, we have that (1 — Zi)Xi-.\ and q i 
are the energy fraction and transverse momentum of the emitted gluon i. We shall use k% and 
qi also to denote the moduli of the transverse momenta (hi = |fc*| and qi = |<fo|). We consider 
the region Zj<l, which gives the leading IR singularity. 

The emission process takes place in the angular ordered region given by 9i > #j_i with 
Qi the angle of the emitted gluon q^ with respect to the incoming gluon k- m . In terms of the 
emitted transverse momenta <& this region is given by 

0i > 9i-i , =► q { > Zi-i^i-i • (6) 

The branching process given in is accurate to leading IR order and, at the inclusive level, 



does not require any collinear approximation ||. The distribution for the emission of gluon i 
is given by 

dVi = — | dzi — A(zi, qi, h) 6 (ft - ^_ift_i) , (7) 

nqf zt 

where angular ordering (|6]) is included. The function A is the form factor which resums 
important virtual corrections for small Z\ 
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In A( Zi , ft, h) = - f dz>m% Q(k t - q') Q(q> - z'q,) . (8) 

J Zi z J q 

This form factor has a simple probabilistic interpretation. It corresponds to the probability 
for having no radiation of gluons with energy fraction x' = z'xi-i in the region X{ < x' < Xj-i, 
and with a transverse momentum q' smaller than the total emitted transverse momentum fcj 
and with an angle 9' > 9, t . Angles and momenta are related by ft — Xi-iE9i and q' ~ x'E9' 
so that angular ordering gives q' > z 'qi. The two boundaries in q' are due to coherence in the 
exchanged gluon (k > q') and in the emitted one (9' > 9i). 

One has 

( 1 k 2 \ 

A(z u qi, hi) = exp [ -a s In — In — ^ , k { > ft, (9) 

so that this form factor has a double logarithmic form, suppressing the radiation both for 
Zi <^1 and for the emitted transverse momentum ft <C kij ' J~zl- 



From the form in (|J) and from the probabilistic interpretation, we have that the function A 
plays a role similar to that of the Sudakov form factor. However there are important differences. 
The Sudakov form factor resums virtual corrections with IR singularities due to soft emitted 
gluons, i.e. powers of ln(l — z), regularising the z — > 1 singularity in the splitting function. This 
implies that when weighted with the energy fraction Zi, the usual branching with the Sudakov 
form factor [Q can be normalised to unity, corresponding to the u = 1 energy sum rule. 

The form factor (§) resums virtual corrections with IR singularities due to soft exchanged 
gluons, i.e. powers of lnx. The branching (|?p cannot be easily normalised to unity for u = 1. 
However this normalisation is not relevant for our study of small x. Note that A(z, q, k) depends 
not only on two transverse momentum scales, as does the Sudakov form factor, but also on 
the energy fraction. This extra dependence in the form factor is one of the important features 
of DIS coherence at small x. A(z, ft k) corresponds in the BFKL equation to the gluon Regge 
form factor, which depends only on k, z and a collinear cutoff (see later). 

Since the intermediate real and virtual transverse momenta are bounded by angular order- 
ing, no collinear cutoff is needed, except on the emission angle of the first gluon. However for 
Zi-i — > one exposes the collinear singularity of ft — > 0. Thus by integrating over the real and 
virtual transverse momenta one generates powers of lnzj_i. After integrating over the energy 
fractions Zj_i one finds that the general perturbative term is of the form 

^ln^x =► -^y, £<n. (10) 



Here each energy fraction integration gives a power of In x while each transverse momentum 
integration gives either lnx or Ink. We have therefore that the contributions with £ > are 
obtained from collinear singularities. Thus if collinear singularities cancel then the leading In x 
contributions are obtained only from IR singularities, i.e. for £ = 0. 

In order to deduce a recurrence relation for the inclusive distribution in the last gluon with 
fixed x = x n and k = k n one has to introduce a transverse momentum p associated with the 
maximum available angle 

n < § => z n q n <p~xE6, (11) 

with xE the energy of the last gluon, k n , which undergoes the hard collision at the scale Q. 
Then one defines the distribution for emitting n initial state gluons 

/n 
J] dVi Q(p - z n q n ) 5(k 2 - k 2 n ) S(x - x n ) . (12) 

The fully inclusive gluon density 

oo 

A(x,k,p) = £ A^(x,k,p), (13) 

ra=0 

satisfies the following recurrence relation 

»2 J — 

A(x,k,p) = A {0 \x,k,p) + /— | — — A(z,q,k) Q{p-zq) A(x/z, \ k + q \,q), (14) 

J itq 1 z z 

where the inhomogeneous term A^°\x,k,p) is the distribution for no gluon emission. This 
equation can be partially diagonalised by introducing the ^-representation 

Au,(k,p) = / dx x w A(x, k,p) . (15) 

Jo 

One finds 



MKP) = A^(k,p) + [^dz^A(z,q,k)Q(p-zq)A UJ (\k + q\,q). (16) 

J 1TQ Z Z 



This equation cannot be further diagonalised in transverse momentum since the kernel depends 
both on the total momentum k and on q and p. Numerical studies are then necessary. 

For the fully inclusive gluon density A(x, k,p) there is a cancellation between the collinear 
singularities which appear in the real and virtual contributions of the kernel. To see this we 
convert the recurrence relation into an inclusive form. By using the identity 

/ 1 ^^^A(^,g,A;)e(p-^) = ^{l - /* dz z» |- A(z, q, k) 6(p - zq)\ , (17) 
JO Z UJ [ Jo oz J 

one finds 

A,(M = A^(k,p) + ^ [*1 [M\k + q\,q) ~ e(k-q)A„(Kqi)}+S u ,(k, P ), (18) 

u J ixq* 



where the inhomogeneous term is 

dq 2 
Q 

and gi = min(q,p). The correction 5uj(k,p) is given by 

d 2 q 



Jg>(k,p) = Ag>{k,p) + H f d ^Q{k-q)Af{hq l ), (19) 

u J q z 



6 u (k,p) = — mM\k + q\,q) 



P -) A(5, ,,*) 



(20) 



where the integration range is q > p. In equation ([181) the first term in the integral comes 
from the real emission contribution while the second is due to the virtual correction. One sees 
explicitly that, as in the BFKL equation, the real and virtual q — > collinear singularities in 
the kernel cancel. Also 5^(k,p) has no collinear singularities since q > p ^ 0. 

For p — > oo the term 8 w {k,p) vanishes and the gluon density A^{k,p) becomes independent 
of p. In fact neglecting 8 w (k,p) and the p dependence in Au(k,p) one finds that ( |18|) becomes 
the BFKL equation for the gluon density T^k) 

72 



?M = j% \k) + 2£ / -| fo(| fc + g |) - e(* - g) *",(*)] . (21) 

uo J Tiq 



Neglecting the p-dependence in Au(k,p) corresponds to neglecting angular ordering. To see 
this we modify the branching distribution in (|7|) and the virtual corrections (|8|) by neglecting 
angular ordering so that the transverse momenta have no lower bound. To avoid singularities 
we have to set a collinear cutoff /i, which, at the inclusive level, becomes irrelevant. The 
modified branching distribution is given by 



nQi Zi (22) 



nqf Zi 

lnA^(z,k) = - j'dz' ^ f^lQ(k-q') Q(q'-fi) 
Jz z J q 



obtained from (|7|) and (|8|) by the substitution Q(qi — Zi-iqi-i) — > ®(qi — /i), and Q(q' — z'q) — > 
Q(q' — p) respectively. This modification has no effect to leading order since the collinear 
singularities cancel. Proceeding as for A, the gluon density satisfies the following recurrence 
relation: 



F{x, k) = J* \x, k) + / ^4 — — A (0) (^, k) Q(q - //) F{x/z, \k + q\), (23) 

J Tcq z z z 



where the inhomogeneous term is related to the one in (|2"T| ) as in (JH5J). From this modified 
branching one obtains the BFKL equation ([21]) in which the momentum has the cutoff p, which 
can be neglected. 



2.1 General properties of gluon distributions 

In the following we discuss some of the properties of the gluon distribution A u (k,p) and its 
comparison with the BFKL distribution !F u (k). As mentioned before, at large p the leading 
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order contribution to Au(k,p) tends to T^ik). We are interested in analysing the subleading 
corrections contained in Au}(k,p) which are due to angular ordering. 

We start by recalling some well known properties of the BFKL distribution. The solution 



^( k ) = L — 13. I U I --77T > ( 24 ) 



of (21) is given by 

dj 1 fk 2 V uf (u;,j) 
2ni k 2 \k 2 ) u - a S x(l) 

where ko and the function /o are fixed by the inhomogeneous term and x{l) is the BFKL 
characteristic function ([3]). 



Taking the initial condition 

„2\ 70 



^<*> = h (I)' 



(25) 



with a given 70 one has that the solution J-^(k) has the same form, and therefore the small-x 
behaviour of T{x, k) is given by 



x -uj c /m 70 
xF(x, k) ~ -j— f p J , w c = a 5 x(7o) • (26) 

For a general initial condition the asymptotic behaviour of T u {k) for k 3> ko and for k ^ k 
is given by the expression (Q) where 7 is the solution of the characteristic equation (|3]) in the 
region < 7 < i and ~ < 7 < 1 respectively. The behaviour of JF(x, fc) at small x is determined 
by the leading singularity of j(as/uj) in the a;-plane which is at 7 C = j(as/uj c ) = 1/2 giving 
u c = aisxih) = 4a 5 ln2. 

We come now to discuss the properties for A(x, k) by taking solutions of the form 

MKp)= ^{^f G{p/kh (27) 

where 7 has to be specified and the function G(p/k) takes into account angular ordering. 
The equation for G(p/k) is obtained by taking the derivative of fll6|) with respect to p: 

pw)=**!, 0(9" ^.«.*) G (rs^n) (t)" • < 28 > 

We consider the case of < 7 < 1 and as a boundary condition we take G(oo) = 1. This 
function depends on as, oj and 7. 

If one takes the initial condition as for the BFKL case, (see (|25|)) 

1 / p\To-i 
^M = ^(rgJ ©(P-*), (29) 

with a given 70, then the solution of the angular ordering equation has the form (p7]) with 
7 = 7q. In this case 7 is a free parameter independent of as and to. 



The expression ( J27Q is a solution of the homogeneous equation, as in the BFKL case, 
provided that 7 is given by the generalised characteristic function which is obtained from (|27| ) 
and (|18D in the limit p — > 00. 



1 = jm,*s) , * = J $ { (^T" o (^ ! - e(* - ,) aim } ■ (30) 

There may be various solutions to this equation and we will consider the leading one, i.e. that 
with the largest uj. In this case 7 is not an independent variable, but is a function of as and 

uj. 

From these equations one finds that the leading order contribution to A(x, k,p) is given by 
J-(x,k). First from ( p8|) one has that the p-dependence of G(p/k) is a subleading correction 
proportional to a$ without 1/uj enhancement. Moreover for G(p/k) — > 1, one has xil-, a s) ~^ 
xirDi the BFKL characteristic function in (|3|). We now list some properties of the angular 
ordering function. 

Behaviour of G(p/k) for p ^> k. In the region q > p ^> k we have A(p/q, q,k) = 1 and 
(fe + q) 2 — q 2 . From (|28|) we have 

"* o &-/*) = *°w/,S(f )"(?)"'■ (3i) 

Since 7 < 1, the derivative vanishes at large p and one finds 

G ^ - 1 - (1 4ffU) iff' 1 • « 32 > 

Behaviour of G(p/k) for p <C k. By using 

In A(p/q,q,k) = -a s \\ia 2 (k/p) Q(k - p) - \n 2 (k/q) Q(k - q) 
and (k + q) 2 ~ k 2 for q ^i k, from ( p8|) one obtains 

^ d 2 q Ik' 
'p nq 2 \q 

where C(u>,as) is a constant which is independent of p. We find then 



pd p G(p/k) ~ a 5 e-^ ln2 ( fc / p ) (£)" |* ^| (-) e ^ ln2 ( fc /«) G{q/k) + C(uj,a s ) 



(33) 



a 5 C(^,a 5 ) fPY „-a s ln\k/p) 



G(p/jfe) ~ J ^ v '^ 7 f e-" s ta ww . (34) 

uj \kj 

This behaviour is similar to that of a Sudakov form factor. The distribution is suppressed 
when the maximum angle 9 available for the initial state emission vanishes, more precisely for 
p ~ xE9 much smaller than the total emitted momentum k. 

Behaviour for 7 — > and as fixed. In this limit the angular ordering function assumes the 
perturbative form, G(p/k) — > Q(p — k). This is as one would expect from fl3"2| ) and (p4|), and 
can be proved directly from (p8f). 

8 



Correction to the characteristic function. The correction 5% (7, as) is given by eq. (|30|) in 
which one substitutes G{p/k) with 5G(p/k) = 1 — G{p/k). From (|3^) one has that 5% = 
xil) ~ X.{li a s0 is regular for 7^0. From 5x one obtains the subleading corrections to the 
gluon anomalous dimension due to angular ordering. By taking into account that 7 is the 
solution of 1 = {ois/u!)x(l) we can write the characteristic equation fl3~0|) in the form 

x(i, a s) = xiri) ■ ( 35 ) 

Expanding ^(7, as) around as = we find 

<9x(7,as) 



7 = 7 + ctsli H , 7i 



X / (7)<9a 5 



(36) 



«c=0 



From these equations we have that the first correction to 7 is as^i ~ a|/a; 2 . To prove 
this, observe that x'(7) ~ 7~ 2 f° r small 7. Moreover, for ct^ — ► and 7 fixed, we have 
Sx(j, cxs) ~ cas where c tends to a constant as 7 — > 0. 

As we shall see from the numerical analysis, the characteristic function ^(7, as) decreases 
with 7, reaches a minimum at 7 C < 1 (for reasonable as), and then rises again. As in the BFKL 
case, we shall denote by Q c the leading singularity in u which corresponds to the minimum of 
the characteristic function at 7 = %. 



3 Numerical methods 



3.1 Evolution in rapidity 

The angular ordered equation is solved by binning the function A(x, k,p) in all three variables. 
To allow the coverage of a wide range of transverse and longitudinal momenta it is convenient to 
store the function on a grid of logarithmic variables y = ln(l/x), In k and hip, where transverse 
momenta are in units of k . This allows us to go to very small x, and to cover the wide range 
of transverse momenta needed to correctly take into account the diffusion in In k. We solve for 
A using the integral equation flUI). From now on we will refer to the gluon density through 
the following function: 



A(y,k,p)=e-yA(e-y,k,p). (37) 



The equation satisfied by A is 



A{y,k,p)=A ( -°){y,k,p) + as / dy' / —±A{e- y ,q,k)Q{lnp + y'-lnq)A{y-y',k + q,q), 

(38) 



where the limits on the q integration, q, \ k + q \ > q min and q,\ k + q\ < g max , are introduced 
because of the finite extent of the grid in Inp and Ink. Our approach is based on the fact 
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that A(0,k,p) = A(°>(0,h,p). Then we attempt to determine A on a grid in y, at points 
y = i5y, where i runs from 1 to some v ax = Y/5y, Y being the highest rapidity in which 
we are interested. Assuming that we know A for all points on the grid up to i5y, then the 
procedure for evaluating the {i + l) th point is the following: as a first approximation we set 



A((i + l)5y,k,p) = A(i5y,k,p). This is put into the integral equation (|38|) to allow us to 
calculate a second approximation to A((i + l)5y,k,p), which can itself be fed in to yield 
a still better approximation. This procedure is repeated until we have a stable value for 
A((i + l)5y,k,p). Generally convergence is reached after about three or four steps. One can 
also aid the convergence by making a better first approximation (e.g. by taking into account 
the first derivative of A with respect to y). In certain regions the form factor A varies very 
rapidly with y', requiring the use of specially adapted integration weights to obtain the correct 
answer. In cases where A varies with y in a known rapid manner, that information can also be 
used. The BFKL equation is solved in a similar manner. 

Later in this article, for the purpose of obtaining precise information about x, it will be 
necessary to go to extremely large rapidities -- y ~ 100, or equivalently x ~ 10~ 50 ! At this 
point the y' integrations, because of their large extent, become very slow, and also require 
one to store information about A at rapidities all the way from to y. This quickly becomes 
prohibitive both in terms of computing time (which scales as y 2 ) and memory requirements 
(one should bear in mind that for each y point, we are storing a large 2-dimensional grid in In k 
and hip, as discussed below). Fortunately the integrand is dominated by small values of y', 
because A(y — y') is a function which decreases exponentially with y' as does the form factor. 
This allows us to truncate the y' integration; a limit of y' < 4:/as is generally found to be 
adequate. 

The other component of the problem is the d 2 q integration. Given that we have stored 
A(y, k,p) on a grid in In k and In p, the task is that of obtaining a discretised kernel K{i k ,i\ k+q \,i q ) 
such that 



'2™ £g,max 



/ ^l A (y,\k + q\,q)= £ £ K(i k ,i lk+ql ,i q )A(y,e^ si ,e^) (39) 

l q — ^9,111111 i \k + q\ 



where k = e lkSe and S£ is the logarithmic spacing between grid points in lnfc and In p. The 
sum over i\k+ q \ is the equivalent of the angular integral. The difficulty that arises is that in 
the region where \k + q\ involves a significant cancellation, a small change in either k or q has 
a large effect on ln|fc + q\ - - so moving slightly away from a grid point defined by i k ,i q , the 
result of the angular integral changes radically. The solution is to think of i k , and i q not as grid 
points, but as extended regions in k and q (and analogously for i\k+q\), and when calculating 
the discretised kernel one must perform an average over the appropriate region. This is found 
to drastically reduce the discretisation errors. 

The main limits on the method described here are due to memory requirements resulting 
from storing the gluon density on a three dimensional grid. Generally the grid resolution 
parameters^ are 5y = 5£ ~ 0.1, together with g min ~ 10~ 6 &o and g max ~ 10 6 A;o. With these 
parameters one can determine A to an accuracy of better than 1%. 



2 Sy and 6£ are kept equal to simplify the treatment of the 6-function in (|38|). 
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As done in section [4.2| , one can impose a dependence of the form k 2 ^~^G^ '(y,p/k) on 
A(°'(y, k,p). The fc 2 ^ -1 ) dependence remains in the solution A, so that it is no longer necessary 
to store the k dimension of A. This allows one to go to smaller bin spacings, increasing 
the accuracy, which is necessary when attempting precision studies of the angular ordered 
characteristic function ^(7, as). To perform the calculation at small 7, one must take into 
account that the integral over the region of small \k + q\ needs to go to extremely small 
\k + q\ <C e -1 /^- 1 (the integral is of the form J dxx 2 ^' 1 ). In principle this would require a 
prohibitively large number of bins - - but the problem has been be resolved by calculating 
analytically the contribution from the region below the lowest stored bin. 



3.2 Iterative method 



We have solved the recurrence relations derived from (|T6|), which in u space read: 

7rg 2 
where the kernel V is given by: 

r w (*,J>, q) = Jo dzz"JA(z, k, q)6(p - zq) 



A(: +1 \k,p)=J^r^k,p, q )A^(\k + q \, q ) , (40) 



a s 



+ A erfc 



+ v^log £| exp [^ + u log \ + a s log 2 | 



(41) 



^ +l \k) = r u (k) [*l0( q -p)J*) (I k + q |) , (42) 

J no 



2^/5? ' V "-0 ~b q2 _ 

with gi = min(p, q), q2 = min(k,p,q) and q^ = min(/c, q). The corresponding equation for the 
BFKL distributions is 

A. 

iiq 2 

where 

r ^) = /W-A W (^) = — 1-Las1 U77 r- (43) 

•>o z w 1 + 2a log ^0(k - //) 

The iterative method is not very efficient (in that it requires more computer memory and 
CPU time than other methods), but has the advantage that it closely mimics the physical 
branching process, thus giving a simple way to evaluate final state (exclusive) quantities, which 
will be examined in a future publication. 

Our method is based on the truncated expansion of the various distributions on a suitable 
basis of Chebyshev polynomials (see also 0). The expansion of the structure function (|40|) 
reads: 

N M 

A^(k,p) « w(k,p) £ £ a& V^[t{k)]Vi M \u{p)\ , (44) 

n=l ra=l 

where w is a weight function, chosen to ensure convergence, the functions t and u map the 
variables k and p (which range from to 00) onto the interval [— 1, 1]. The efficiency of the 
method relies heavily on a suitable choice of the function t and u. We have used the usual 
logarithmic mapping: 

/-(/,-) = I ln\/^ I In k , u(p) = fln,/^) In V (45) 
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where k min , p min and /c max , Pmax are respectively the lower and upper cutoffs in the k and p 
variables. 



The basis functions are defined by 

nv-i 



vl N \t) = | 



y: T 4 (t n )T,(t) + ~ 



(46) 



where Tj are Chebyshev polynomials, and are an orthonormal complete set on the points {t n }, 

tn = C ° S ( ~2N + n ) ' n = 1 '---' iV ( 47 ) 

with the notable property that V^ N \t m ) = 5 n , m - 

With these definitions, the expansion coefficients are just the (rescaled) function values on 
the N x M rectangular grid {(k n ,p m )}: 



(r) __ A^(k n ,p 
w{k n ,p 



m ) 
m ) 



(48) 
where t(k n ) = t n , n = 1, . . . ,N and u(p m ) = u m , m — 1, . . . , M. 



Inserting expansion (44) , the recursion ( ^0|) becomes 



N M 

a, 



(r+1) - E E WwaJL, (49) 



nm 

n'=l m'=l 



where the matrix 



L nm y m > = f — T U ;i + q [ q) ^(k n ,p m , q)V^[t(\ K + q \)}V^\u{q)} (50) 

must be evaluated only once, and then used to iterate the equation as many times as desired. 

We have checked the stability and the convergence of this method by varying the functions 
w, t and u, and the parameters N and M. The method is weakly sensitive to the choice of w, 
as long as the starting condition ,^ 'f is smooth, and numerical convergence is ensured for 
the integral appearing in ( |50|) . We have used w(k,p) = A^°'(k,p) for the Aw distribution and 



7T 



w(k) = J jj\ 2 fc> r the J-'w case. The choice ([45]) allows us to obtain good stability and fast 
convergence with iV and M of the order of 30 — 40. 

The number of iterations needed to obtain a stable solution depends drastically on the value 
of uj. Far from the u;-plane singularity (u ~ 1), 20 — 30 iterations are sufficient to obtain an 
accurate solution over all the /c-range (the iteration converges faster for smaller values of k). 
On the other hand, the required number of iterations increases dramatically as u approaches 
the critical value u c : with uj — uj c w 0.01 we need about 5000 iterations to obtain a reliable 
result. 
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Figure 2: Plot of A(x, k,p)/A(x, k,p) against k for different values of p; x = e 10 , p = e l5 k 
and as = 0.2. 



4 Numerical results 



In this section we report the results obtained by solving ([14]) and (|T6|) for the gluon distributions 
A(x,k,p) and A^ik.p) respectively. We compare these results with those obtained from the 
BFKL equation. 



4.1 Behaviour at small x 



First we study A(x, k,p) for a simple initial condition 

A i0 \x, k,p) = 5(1 - x) \8{k - k Q )Q(p - q v 

k 



(51) 



where k sets the momentum scale. As a collinear cutoff for the first emission we take ^mm- 
This condition is not quite physical but is suitable for studying the general properties of the 
solution. We first show that A(x, k,p) becomes independent of p for p 3> k. Then we show 
that its behaviour for x — > is less singular than that of the BFKL gluon distribution. 

p- dependence. In fig. 2 we plot A(x, k,p) as a function of k for increasing values of p and for 
fixed x and as- As p increases the gluon distribution becomes independent of p. To show this 
we plot the ratio A(x, k,p)/A(x, k,p) with p in the asymptotic region. As expected from the 
discussion in sect. 2 the limiting value is obtained first at low values of k. 

uj-plane singularity. In fig. 3 we plot the two distributions A u {k,p) and J- U (k) as a function of 
uj at fixed k = ko and as- The value of p in Au}(k,p) is in the asymptotic region. From (^) one 
has that !Fu(k) diverges at the singular point u c = a,sx(l/2) = 4 In 2ag, the minimum of the 
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Figure 3: A U} (k,p) and J- U {k) for as = 0.2 and p = e 15 k . For the BFKL case the dotted line 
corresponds to the exact singularity at u = u c = 0.5295. For the CCFM case the dotted line 



corresponds to a singularity at to = tu c = 0.4301, as determined in section W72 



BFKL characteristic function at 7 = 1/2. The gluon distribution Au(k,p) has a singularity at 
a value[] Cj c ~ 0.4301 smaller than uj c ~ 0.5295, the BFKL value. 

For the BFKL case the numerical distribution actually tends to diverge at a value of u 
which is 2 — 3% smaller than u c . This is due to the presence of lower and upper limits, g m i n 
and g max , in the q and k grid used for the numerical calculation. It is known || that if the 
transverse momentum range is finite, the ^-singularity is shifted to a lower value of the order 
ijj c — f u c (l — vr 2 /ln 2 g min /g max ). Similar behaviour has also been noted |TIJ in the context of 
the dipole approach to small- x physics JTT[] . For the values of g min and g max considered here 
one finds 7r 2 /ln q m m/q max . — 0.01. The reason for this shift is that BFKL diffusion from the 
edges of the grid modifies the shape of the fc-distribution which in turn leads to a reduction 
in the observed power. With angular ordering, the diffusion is reduced (this will be discussed 



also in section 4.2) and therefore the edges of the grid have a smaller effect. 



Small-x behaviour. Another way to obtain the position of the cu-plane singularity consists in 
studying the small-x behaviour of A(x, k,p) and T{x, k). In particular in fig. 4 we plot in the 

small-x region 

d 11 

— lnx.F(x, k) ~cu c - — , y = ln-, (52) 

dy 2y 



x 



where —\j2y is the first subasymptotic contribution. One expects that there might be a similar 
kind of behaviour in the case of A(x, k,p), with uo c replaced by uo c . This is indeed seen and 
as before we find u c < uj c . The values for u c and uj c agree with those found in the previous 
analysis in fig. 3. In the BFKL case we have plotted the analytical result obtained for g min — > 



3 The value quoted here is actually the one determined in section 4.2, which has a higher precision than that 
obtained by examining the position of the singularity of A^ . 
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Figure 4: The effective power as a function of 1/y; as = 0.2 and p = p = 3 x 10 e 



and g max — > oo. With finite values of g min and q max , increasing y, the width in k of the solution 
increases until it is comparable to the extent of the finite grid, at which point the numerical 
curve flattens off. This is the same phenomenon that was noted previously for the shift of the 
singularity of T&. For the angular ordering case, we have plotted an analytical line analogous 
to the BFKL one, with the asymptotic power £u c = 0.43 fitted to give agreement with the 
numerical results. The tailing off of the numerical curve at large y is also observed in the 
angular ordering case, but it sets in later than for the BFKL curve — consistent with the idea 
that diffusion is reduced by angular ordering so that the width of the solution approaches that 
of the grid only at larger y. 

We have also studied the /c-dependence for large y. In the BFKL case the numerical result 
fits well the expected behaviour ~ l/kk , corresponding to 7 = 7 C = 1/2. With angular 
ordering we find A(x,k,p) ~ k 2 ^^ 1 ^ with 7 = 7 C ~ 0.61. We shall analyse this behaviour in 
more detail later. 



4.2 Characteristic function 



We now study the generalised characteristic function \ an d the corresponding angular ordering 
function G(p/k). Recall that for small x the gluon distribution A(x, k,p) has the asymptotic 
form 

xA(x, k,p) ~ - 1 - r ( -^ G{p/k) , uj c = a s xil, as) . (53) 



k 2 



n'n 
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To obtain x an d the corresponding angular ordering function G(p/k) as a function of 7 and 
a 5 we use the following method. We solve equation (JTJJ) by using a trial initial condition 



A (0 \x,k,p) = -( T1 ) 6(l-x)Q(p-k), (54) 




with a fixed value of 7 and as- From the discussion in section 2 one has that asymptotically 
for x — > 0, A(x, k,p) has the form (|53|) with the same 7 as the inhomogeneous term ([54]). The 



x- and p-dependence of the initial condition is not important. Since the solution has the form 
of ( |53| ) we only need to deal with G, which depends just on p/k. This means that one doesn't 
need to store the k dependence of the solution and the finite extent of the grid no longer has 
a significant effect, drastically reducing the errors. 

Characteristic function x- To obtain ^(7, as) by using the initial condition (|54[) we compute 
A(x, k,p) with p in the asymptotic region. By taking the small-x limit we determine the 
intercept uj c with a high accuracy — the relative error is of the order of 10 -4 for much of the 
7 region. From u) c we obtain the characteristic function 

X(7.as) = — , (55) 

a s 

as a function of 7 and as- In fig. 5a we plot x as a function of 7 for various as- We plot for 
comparison the BFKL characteristic function x- We see that 6x = X ~ X is positive, increases 
with 7, and increases with as- Moreover we find 5x ~ 7 for 7 — » (as small and fixed) and 
5x ~ «5 for as — > (7 small and fixed). This agrees with our earlier observation in section 2 
and implies that the next-to-leading correction to the gluon anomalous dimension coming from 
angular ordering is of order a|/o; 2 . In fig. 5b we plot Sx/c^s which shows that there are notable 
next-to-next-to leading corrections especially at large 7. 

The symmetry for 7 — > 1 — 7 of the BFKL characteristic function is not valid for x- 
Recall that this symmetry is based on the fact that the regions of small and large k are equally 
important. In the CCFM case, angular ordering favours instead the region of larger k. However, 
small values are still accessible. Therefore the function x decreases faster than x for increasing 
7, but, after a minimum at a point 7 C larger than the BFKL value 1/2, x increases again. In 
fig. 6a and 6b we plot as a function of as the values Xc and 7 C with Xc the minimum of x an d 7 C 
its position. As expected the differences compared to the BFKL values Xc — 4 In 2 and 7 C = 1/2 
are of order 0:5. These results are consistent with the asymptotic solution in fig. 4. From our 
determination here, we obtain a very accurate estimate of the position of the singularity in 
u: for as = 0.2, we find u) c ~ 0.4301 (the corresponding BFKL value is uj c ~ 0.5295) and 
7 C ~ 0.6106. 

Figure 6c shows the second derivative, x"i of the characteristic function at its minimum; this 
quantity is important phenomenologically because the diffusion in In k is inversely proportional 
to the square root of x"- From the figure, one can see therefore that the inclusion of angular 
ordering significantly reduces the diffusion compared to the BFKL case. 



Angular ordering function. By using the trial initial condition (|54[) and taking the small-x limit 
for A(x, k,p) we compute the angular ordering function G(p/k) at the given value of 7 and as 
(see eq. (|53D). In fig. 7 we plot G(p/k). The behaviours of G(p/k) for small and large p (see 
(|32| ) and (|34|)) are as expected. In particular, we note that as as — > (fig. 7a) G{p/k) tends 
slowly to become 1 everywhere; as 7 — ► 0, G(p/k) tends towards the function Q(p — k). 
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Figure 5: (a) The characteristic functions with and without angular ordering; x{li as) and 
x{l) are plotted as functions of 7 and 7 respectively, (b) The difference, 5x = x{l) ~ x{li a s), 
between the BFKL and angular ordered characteristic functions, divided by as- 
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Figure 6: (a) The value of the minimum of the characteristic function, Xc as a function of as- 
(b) The position of the minimum of the characteristic function, %, as a function of as- (c) 
The second derivative of the characteristic function, x", at its minimum, as a function of as- 
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Figure 7: The angular ordering function G(p/k): (a) for a range of as, (b) for a range of 7. 
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5 Conclusions 



In this paper we have studied the contributions to the subleading corrections of the small- x 
gluon density which are due to angular ordering. Since they are based on the physical property 
of QCD coherence, one expects that they are among the important corrections. Another 
important subleading contribution, that which fixes the scale of the running coupling ||, is 
not included in our study. The calculation has been done mostly by numerical methods, which 
prove to be quite reliable. In future papers they will be extended to the study of associated 



distributions |12|, [13], [14] for which angular ordering is relevant already to leading level. 



Our main results are summarised in figs. 5-7 in which we plot the generalised characteristic 
function xili a s) and the angular ordering function G{p/k). From these plots we have studied 
the subleading corrections 5%(7, ots) = x{l)~ x{li a s) and 5G{p/k) = G(oo) — G(p/k). We find 
that xili a s) decreases with 7 faster than the BFKL characteristic function, it has a minimum 
at 7 C which is larger than 7 C = 1/2, the BFKL critical point, and it rises again at larger 7. 

The angular ordering function G(p/k) has the structure of a typical form factor: it vanishes 
when the maximum available angle 9 vanishes, i.e. for p ~ xEO much smaller than k, the total 
emitted momentum. 

The BFKL symmetry 7 —>■ 1 — 7 is lost since conformal invariance is broken by angular 
ordering. The physical basis of conformal invariance is that in the BFKL equation the regions 
of small and large momentum are equally important. The coherent branching instead tends to 
evolve toward large momenta. However, at each intermediate branching, the region of vanishing 
momentum is still reachable for x — > 0. Within the angular ordering formulation, this effect 
has been discussed also in W% . 



The fact that during the branching the intermediate momentum could vanish implies that 
the evolution contains non-perturbative components in an intrinsic way, not only in an initial 
boundary condition. It should be noted that in this non-perturbative region the distribution is 
non-singular (collinear singularities cancel), so that for any small but finite x non-perturbative 
effects of the small-A; region are not too important. However they become very important for 
the asymptotic limit, x — > 0. As in the BFKL case the small- fc region generates a singularity 
in u at Cb c > and an anomalous dimension 7 C which is non- vanishing with as (see fig. 6b). 

By expanding ^(7, as) in powers of as at fixed as/uj one obtains the part due to angular 
ordering of the next-to-leading correction asli^O-s/^) (see (0)). We have not obtained its 
analytical form but we have shown that the angular ordering correction in the small- 2 limit 
starts with a power not smaller than a s /u 2 . 
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